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C| Abstract. We show that electron tunneling from edge states in two-dimensional 

| topological insulator into a parallel electron waveguide leads to the appearance of 

O ■ spin-polarized current in the waveguide. The spin polarization P can be very close 

to unity and the electron current passing through the tunnel contact splits in the 
waveguide into two branches flowing from the contact. The polarization essentially 
depends on the electron scattering by the contact and the electron-electron interaction 
in the one-dimensional edge states. The electron-electron interaction is treated within 
'"^J ■ the Luttinger liquid model. The main effect of the interaction stems from the 

. renormalization of the electron velocity, due to which the polarization increases with 

the interaction strength. Electron scattering by the contact leads to a decrease in P. 
A specific effect occurs when the bottom of the subbands in the waveguide crosses 
the Dirac point of the spectrum of edge states when changing the voltage or chemical 
■ potential. This leads to changing the direction of the spin current. 

r< 

00 
(N 

OV 

o: 

(N 



X 



o 
o 



PACS numbers: 72.25.Dc, 73.40. Gk, 73.63.Nm 



Spin current in an electron waveguide tunnel- coupled to topological insulator 2 
1. Introduction 

Spin properties of surface and edge states are hallmark of topological insulators 
(TIs) [H |2] which attract much interest as a playground for fundamental physics and 
a basis for potential applications in spintronics and quantum computing. The surface 
states in three-dimensional TIs and edge states in two-dimensional (2D) TIs have a 
helically polarized spin structure. The spin of electrons is strongly coupled to their 
momentum. A consequence of this coupling is that the edge states carry a spin current. 
This fact stimulates the idea of creating spin currents in a nonmagnetic conductor, 
which is in contact with a TI. Indeed, if the tunneling electrons conserve their spin 
and direction of the velocity, the spin filtering of electrons occurs in the nonmagnetic 
conductor since the electrons with different spins move with different velocities. 

The existence of the helical edge states with Dirac dispersion in 2D TIs was proved 
by recent experiments [31 HJ [5J [6]. Rapid progress in nanotechnology allows one to 
expect the creation of mesoscopic heterostructures containing TIs [SIC]- Therefore it is 
interesting to consider a model system in which a 2D TI is coupled to a strip of normal 
2D electron gas via a tunnel contact. Tunneling of electrons from the edge states in 2D 
TI into the normal conductor results in the transition of electrons with different spins 
in the states with different momenta. In that way a spin current arises in the normal 
conductor even if the barrier is nonmagnetic and external magnetic field is absent. 

The present paper aims to study the spin polarization of the tunnel current in 
the normal conductor, to clarify what factors affect the polarization and estimate the 
polarization for realistic conditions. The edge states are one-dimensional (ID) and 
therefore electron-electron interaction plays an essential role in their tunneling. The 
ID electrons in these states can be considered as a helical Luttinger liquid [5]. The 
interaction is known to change the one-particle density of states and to renormalize 
the electron velocity [9]. Both these factors affect the tunnel current. The tunneling 
between parallel ID systems is a subject of high activity in last decade. Semiconductor 
heterostructures with tunnel-coupled quantum wires were created, and impressive 
magneto-tunneling experiments were carried out |10j . which have revealed many-particle 
effects. The theoretical works, which are closest to the present paper, have studied 
the electron-electron interaction effects on the tunnel current [TT], the dependence of 
the current on the contact length [12], and such many-particle effects as fractional 
quasiparticles and spin-charge separation [13]. Spin transport has been investigated 
much less. Spin filtering of the tunnel current was explored in the case where the spin- 
orbit interaction acts in one of parallel quantum wires [H]. Recently the theory of 
tunneling between ID systems was extended to the helical Luttinger liquids [IS]. Spin 
injection from helical edge states into a normal electron gas was not studied. 

We solve this problem by considering a model system consisting of a 2D TI coupled 
to a strip of normal 2D electron gas via a tunnel contact of finite length. Electrons 
in the edge states are treated as helical Luttinger liquid and the strip is considered 
as a multimode electron waveguide. We have found that the electron current flowing 
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through the tunnel contact into the normal conductor is split there into two currents of 
electrons with opposite spins flowing in opposite directions. Thus, the electron current 
in each branch of the waveguide outside the contact equals half of the total current 
through the tunnel contact. It is interesting that the electron current in the normal 
conductor is highly spin polarized and there is a spin current. The spin polarization 
P can be very close to unity. The polarization essentially depends on the interelectron 
interaction in the Luttinger liquid and the electron scattering by the contact. The main 
effect of electron-electron interaction is caused by the renormalization of the electron 
velocity in the edge states of the TI. Increasing the interaction results in the growth of 
the polarization. The scattering by the contact leads to a decrease of P. A specific effect 
arises when the applied voltage is varied so that a subband bottom in the waveguide 
crosses the Dirac point of the edge state spectrum. Under this condition the sign of the 
spin current is changed. 

The paper is organized as follows. Section [2] describes the model we use for 
calculations and the qualitative discussion of the main issues. The formalism of 
calculations is presented in Sec. |3j Main results are presented and discussed in Sec. HJ 
Sec. [5] contains concluding remarks. 

2. The model 

To study the spin-polarized current, which can be created due to the tunneling of 
electrons from a 2D TI, we consider the model system shown in Fig. [TJ The TI occupies 
half-plane y < 0. At the line y = 0, it is coupled to a 2D strip of normal electron gas via 
the tunnel barrier of length L. A voltage V is applied between the TI and the strip so 
that the chemical potential of the edge states fijj is higher than the chemical potential 
in the strip /j,^. Electrons in edge states have gapless Dirac-like spectrum, with the spin 
being strongly coupled to the momentum p: the spin of the right-moving electrons is 
directed upwards and the electrons moving to the left have a spin down. In the strip 
of width W, the electron spectrum is quantized in transverse direction. Within each 
subband the spectrum is parabolic in the momentum k. The choice of a normal electron 
gas, i.e. an electron gas without spin-orbit interaction, as a medium into which electrons 
tunnel from the TI, is caused by the fact that here the spin current is uniquely defined 
and can be measured experimentally (see a discussion in Ref. [16] and references there). 

Spin current depends on the energy of tunneling electrons. Qualitatively the 
tunneling process looks as follows. At the beginning let us neglect the scattering by 
the contact. If the energy is larger than the Dirac point of the edge state spectrum, 
right-moving electrons with spin up pass from the TI into right-moving states in the 
strip. In contrast, the left-moving electrons with spin down come to the left-moving 
states. Thus the tunnel current is split in the strip into two opposite electron flows with 
opposite spins, as it is shown in Fig. [TJi,b. In that way the spin current arises. 

If the energy is below the Dirac point, the picture is changed as shown in Fig. [lb. 
The right-moving electrons in the TI tunnel into the left-moving states in the strip 
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Figure 1. (Color online) (a) Tunnel contact of 2D TI and the strip of 2D normal 
electron gas. (b) Energy spectra of edge states (straight lines) in the TI and electrons 
in the strip (parabolas) in the case where the first subband is above the Dirac point 
and [Iti > (In > 0. Occupied states are shown by heavy lines, (c) Energy spectra in 
the case where the first subband is below the Dirac point and fiN < fJ-Ti < 0. 



and vise-versa, the left-moving electrons in the TI come to the right-moving states. 
Correspondingly, in this case the spin current has opposite direction. 

Thus, the shift of the chemical potentials of the TI and the normal electron gas, fj,jj 
and /An, through the Dirac point results in the change of the sign of the spin current. 

In a realistic situation, there are other factors, which strongly affect the polarization 
in addition to the conservation of the spin and velocity direction discussed above. 
These are the scattering of electrons by the contact, due to which some electrons can 
change the direction of their velocity [HI [15], and the electron-electron interaction, 
since it substantially changes the spectral function of electrons in ID edge states and 
renormalizes the electron velocity [9| [15]. 

3. Formalism 

The total Hamiltonian of the system is 

H = H E s + H N + H tunn , (1) 

where Hes is the Hamiltonian of the isolated edge states, the Hamiltonian of the 
isolated strip, and H tunn describes the tunneling between the edge states and the strip. 
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Electrons in the edge states can be considered as a helical Luttinger liquid [S] and 
described by bosonized Hamiltonian [T7] 



Hes = -y / <!■>■ 



it 8 (x) + — 



g 2 \ dx } ' ' 

where and H(x) are conjugated bosonic phases, vf the Fermi velocity, and g the 
interaction parameter. The ID fermion creation operator is presented in terms of chiral 
fields ip ESr for right and left-moving fermions (r = ±1) 

4;s = X^W> ( 3 ) 



r=± 



with ip ESr being 



, f _ / irppx 

' ES * ~ V 271V 6 



exp |i\/7r r6>(x) + ^ cfa/n(a/) | 



(4) 



Here v = Vp/g is the electron velocity renormalized by the interaction, pp = nu/vp, 
ftw c is the energy cutoff. 

Electrons in the strip are supposed to be non-interacting gas described by the 
Hamiltonian 

Hn = Z2( ~2^" + £ n + E cj 4,n«fc,n > ( 5 ) 

where e n is the quantization energy of nth subband, a\ n and a^n are creation and 
annihilation operators, E c is the energy of the conduction band bottom in the infinite 
2D electron gas. 

The tunneling Hamiltonian is conveniently written in momentum presentation, 
i.e. in terms of the operators a\ n , ak, n and Fourier harmonics of the field operators 

+ H.c. , (6) 

p,k,n 

with c)p being the creation operator for ID fermion with momentum p. The tunneling 
matrix element is easy to find in the case where the strip is a square potential well of 
the width larger than the penetration length of electrons under the barrier, 

it (v k)? ~ r 2 M2 / wfrzjwa y (7) 

\T n {pM ~T lNhsW3 [ (p _ fc)i/2 ) > W 

with In and lps being the normalization lengths for the strip and edge states. 

The electron and spin currents in the system are compactly expressed in terms of 
a matrix I r ' 1 defined as a tunnel flow of electrons from chiral state (r = ±) in the TI 
to the right- and left-moving states (I = ±) in the strip. We calculate I r ' 1 in the first 
approximation in H tunn . The standard procedure [18] gives an expression for J r,/ in 
terms of spectral functions and the Fermi distribution function /(e) = 1/[1 + exp(e/T)]: 

^= { EE\Tn(p,k)\ 2 f£A n (k,e+eV) 

lk>0p,n (8) 

xA ES ( P ,e)[f(e)-f(e+eV)). 
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Here A n (k,e) is spectral functions of non-interacting electrons in the strip 



An(k,e) = 2-kS e 



£n — E c + fi TI ) ; 



h 2 k 2 

2m 

A r ES (p, e) is the spectral function of the chiral ID fermions [11 



e—rhv(p—rpp) 



e+rhv(p— rpp) 



+(e — > — e, p — rpi? — > —p + rpp) 

where a = (g + 1/g — 2)/4, 



(9) 



(10) 



1 2 + 2^T 



11^ 



fej c r(a) V2ttT 
with r(a) being the gamma function. 

The tunnel current J through the contact is given by 

r=± l=± 

The electron currents J 1 in the right (/ = +1) and left (I = —1) branches of the strip 
outside the contact are determined as follows 



j 



r =± 



(13) 



(14) 



The spin current J l s in the right and left branches of the strip reads 

The spin polarization of the current in the strip is defined as the ratio of the spin over 
the charge current in each branch 



p _2e4 



(15) 



From the symmetry of the matrix I r ' [see Eq. OH])] it is evident that P^i = —Pi. 



4. Results and discussion 



We have calculated the electron and spin currents and studied in detail their dependence 
on such parameters as the voltage drop on the tunnel contact, the size of the contact, 
the interaction parameter in the edge states and the chemical potential of the system. 

The electron current J through the contact as a function of the applied voltage V 
is shown in Fig. [2] for small voltage (eV < firi)- In what follows, the dimensionless 
notations are used: the current, the spin current, the length and the energy are 
normalized correspondingly to J = 2em 3 VpT 2 L/7ch 5 , J s = Jh/2e, I = h/2mvp, 
e = 2mvp. The voltage is normalized to e/e. 

The results presented in Fig. [2] were obtained in the case of long contact Lp F ^> 1 
(Pf being the Fermi wavevector in the edge state) where the scattering of the momentum 
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Figure 2. (Color online) Tunnel current through the contact as a function of the 
voltage for different strip width W and the interaction parameter g. (a) W = 301, 
g = 0.9 (line 1), g = 0.6 (line 2). (b) W = 5001, g = 0.9 (line 1), g = 0.6 (line 2). 
Other parameters are: E c = 0, fi = 0.5s, T — O.Ole, L = 3001. 



is small and therefore the polarization is close to unity. The current oscillations in Fig. [2^ 
are caused by the quantization on the strip width. Increasing the interaction strength 
shifts the first peak to higher voltage. The peaks of the current are smoothed out with 
increasing temperature and decreasing g and L. 

In the case of the wide strip, W 3> W c = AttH/ ■sjm ■ max(eV, T), the oscillations 
disappear (Fig. [2b) and the current- voltage characteristic turns into a power dependence 
(Fig. [2b) that is standard for tunneling in the Luttinger liquid [9] . 

A similar oscillatory behavior occurs also in the dependence of the current on the 
chemical potential. 

The most interesting result is obtained in calculations of the spin polarization of 
the current in the strip. Fig. [3] shows the spin and electron currents as a function of 
voltage for a rather wide strip (W 3> W c ). At small voltage, the spin current practically 
coincides with the electron current, however with increasing voltage the spin current 
drops to zero and then changes its sign. Correspondingly, the polarization changes 
from unity to a negative value, which goes asymptotically to -1. This behavior of P 
is explained as follows. When the voltage is small, the energy layer eV, within which 
current- carrying electrons tunnel from the TI into the strip, lies above the Dirac point. In 
this case, practically all electrons with spin down come into left-moving states and spin- 
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Figure 3. (Color online) Spin polarization, the charge and spin currents in the left 
branch of the strip as a function of the voltage. The used parameters are: W = 2001, 
L = 500f, E c = 0, uti = O.le, T = O.Ole, g = 0.9. 




Figure 4. (Color online) The dependence of the polarization on the electron-electron 
interaction parameter for various voltages. The used parameters are: W = 2001, 
L = 500/, E c = 0, uti = O.le, T = O.Ole, g = 0.9. 



up electrons come to right-moving states in the strip (see Fig. [T]d) and therefore \P\ ~ 1. 
With increasing voltage, the above layer of tunneling electrons expands downward. 
When its lower boundary intersects the Dirac point, a fraction of the tunneling electrons 
with spin up comes into the left-moving states (see Fig. [Tfc) and similarly some spin- 
down electrons come to the right-moving states. As a result, the polarization decreases 
and changes the sign. 

A similar effect takes place also in the case of the polarization dependence on the 
chemical potential. The polarization changes the sign when \xti is decreased below the 
Dirac point. The width of the transition interval increases with increasing V and T and 
with decreasing contact length L. 

An interesting question is how the electron-electron interaction strength affects the 
polarization. The interaction effect is described by the parameter g, which decreases 
from 1 to on increasing the interaction potential from zero to infinity. Within the 
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Figure 5. (Color online) The electron and spin currents as a function of the contact 
length, W = 150Z, fi T i = 0.12e, E c = 0, T = O.Ole, g = 0.9, eV = 0.05e. Inset: the 
dependence of the polarization on the interaction parameter for various lengths L. 

Luttinger liquid model, the decrease in g has two consequences that could be important 
for the current polarization. First, the velocity of the chiral fermions increases and 
second, the ID fermions are distributed over wider region of momentum space. 

The dependence of P on g is presented in Fig. HJ The polarization is seen to increase 
with decreasing g. This fact is explained by the first of above factors, namely, by an 
increase in the velocity. Actually, increasing the velocity leads to a shift of the Dirac 
point down and therefore a fewer electrons fall into the energy interval below the Dirac 
point. Fig. H] shows also that at high enough voltage the polarization sign changes with 
increasing g. 

Now consider the dependence of the polarization on the length of the contact, 
which reflects the effect of the electron backscattering by the contact of small length. 
Fig. [5] shows the dependence of the charge and spin currents in the strip on the contact 
length. As the contact length is small (Lpp <C 1) the charge current grows with L 
as J ~ L 2 [19] and the spin current is negligibly small because of strong scattering 
by the contact. When Lpp > 1 the spin current increases, because the scattering 
becomes weaker, and asymptotically reaches the electron current at Lkp ^> 1, with 
kp = ^j2m{^Ti — E c — £i)/h. In this limit, the currents J and Js increase linearly with 
L until the perturbation theory in the tunnel Hamiltonian is not disturbed [13]. The 
inset shows that the polarization decreases with decreasing interaction strength similarly 
to Fig. H 

5. Conclusion 

We have studied the electron tunneling from helical edge states of a 2D TI to a normal 2D 
electron gas and concluded that this process can be an effective mechanism of creating 
the spin current in normal conductors. 

The mechanism consists in the fact that the electron tunneling into the normal 
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conductor gives rise to two oppositely directed currents of electrons with opposite spins. 
One current is directed to the right from the contact, and the other does to the left, so 
that an single continuous spin current flows parallel to the common boundary of the TI 
and the normal electron gas. 

We have studied in detail this mechanism in detail for the system where a 2D strip 
of normal electron gas is coupled to the 2D TI via a tunnel contact of finite length. The 
electron current along the strip is equal half of the tunnel current through the contact 
and is highly polarized. The spin polarization essentially depends on the contact length, 
the electron-electron interaction parameter of the edge states and the voltage drop across 
the contact. 

The most significant features of the spin-polarization behavior are as follows: 

(i) the polarization increases with increasing contact length on the scales Lpp ~ 1 or 
Lkp ~ 1; 

(ii) the direction of the spin current changes with changing voltage or the chemical 
potential when the bottom of the subbands in the strip crosses the Dirac point of 
the edge state spectrum; 

(iii) the polarization grows with increasing electron-electron interaction strength due to 
the renormalization of the electron velocity caused by the interaction. 

The effect we have described here could be observed experimentally, for example, by 
measuring the spin accumulation near contacts to the strip. So, if the strip is connected 
with electron reservoirs containing localized spin states, which are a spin bath, the spin 
density is accumulated there, in the same way as has been shown recently in Ref. |20j . 
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